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Objectives 

Almost  every  aspect  of  Command  and  Control  (C2)  deals  with  making  decisions  in  the 
presence  of  uncertainty.  The  uncertainty  may  come  from  noisy  data  or,  indeed,  regions  of  the 
battlespace  where  there  are  no  data  at  all.  Statistical  models  that  accoimt  for  noisy  data  are  well 
accepted  by  the  science  and  engineering  commimity,  but  the  full  quantification  of  uncertainty 
due  to  lack  of  knowledge  (caused  either  by  hidden  processes  or  missing  data)  is  a  challenging 
statistical-modeling  problem.  When  the  hidden  process  is  a  random  map,  the  problem  is  even 
more  challenging. 

A  map  is  a  powerful  way  of  summarizing  spatial  data  and  information,  and  our  view  is  that 
from  great  maps  come  knowledgable  command  decisions.  Different  commanders  have  different 
“apertures”  and  hence  the  maps  should  depict  information  consistently,  regardless  of  scale.  The 
battle  commander  needs  more  global,  aggregated  information,  whereas  the  platoon  commander 
is  often  making  decisions  based  on  local,  more  focused  information.  Our  emphasis  is  on  making 
maps  that  are  statistically  optimal  and  are  easily  interpretable  at  different  resolutions/apertures. 

Approach 

We  took  two  approaches.  The  first  approach  combines  the  two  aspects  of  noisy  or  missing 
data  and  change-of-aperture,  to  construct  statistical  models  that  deal  directly  with  computation  of 
posterior  statistical  distributions  of  the  map.  Because  of  the  hierarchical  nature  of  C2  problems 
involving  change-of-aperture,  models  can  often  be  built  using  an  acyclic  directed  graph  (ADG). 
In  an  ADG,  each  ‘node’  can  have  ‘parents’  that  it  depends  on,  and  the  node  can  have  ‘children’ 
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that  depend  on  it.  If  each  node  in  an  ADG  has  only  one  parent,  then  the  ADG  is  a  tree.  Spatial 
dependence  can  be  incorporated  through  the  ADG’s  parent-child  relationship,  which  is  a  type  of 
multi-resolution  model.  In  the  case  of  a  tree-structured  model,  very  fast  statistical  algorithms  can 
be  developed  for  removing  noise  and  filling  in  missing  data  at  all  apertures.  The  resulting 
posterior  distributions  are  summarized  in  different  ways,  depending  on  the  command  queries. 

In  the  second  approach,  we  developed  statistical  methodology  to  produce  an 
(approximately)  optimal  map  for  all  queries.  Most  queries  are  nonlinear  in  nature,  but  the  most 
fi'equently  used  optimality  for  map  estimation  is  geared  toward  answering  linear  questions.  The 
approach  we  take  is  to  produce  a  single  map  estimate  that  preserves  certain  first-  and  second- 
order  moment  properties  of  the  underlying  model,  making  it  approximately  optimal  for 
answering  nonlinear  questions.  One  particular  map  we  concentrated  on  was  the  danger  field. 

Finally,  the  temporal  component  has  been  incorporated  to  allow  fast  updating  of  spatial 
maps  based  on  newly  acquired  C2  data.  As  a  result,  the  maps  can  be  animated  and  their 
evolution  studied  during  the  course  of  a  battle. 

Results 

The  research  focused  on  the  mapping  of  the  danger  field  and  like  quantities  in  the 
battlespace.  A  hierarchical  framework,  along  with  an  object-oriented  computer  program,  has 
been  developed  to  describe  the  spatio-temporal  evolution,  and  interaction,  of  battlespace 
constituents  (Cressie  et  al.,  2002).  The  program  simulates  noisy  data  coming  from  different 
sources  and  at  different  spatial  and  temporal  resolutions.  The  data  are  used  to  map  the  danger 
potential  of  the  enemy  at  any  location  in  the  battlespace.  The  resulting  map  can  then  be  used  to 
answer  diverse  queries  at  different  levels  in  the  chain  of  command.  Due  to  the  diversity  of  the 
queries,  care  needs  to  be  taken  in  obtaining  a  statistically  optimal  map.  The  approach  taken  in 
Aldworth  and  Cressie  (2002),  covariance-matching  kriging  (CM),  is  intuitively  appealing  and 
computationally  fast.  Cressie  and  Johannesson  (2001)  give  a  practical  implementation  of  the 
CM  method  in  spatial  settings.  Also,  a  paper  is  in  preparation  giving  a  fully  Bayesian  approach 
to  estimating  the  danger-potential  fields  using  sequential  imputation. 

To  deal  with  consistency  of  predictions  across  changes  of  aperture,  we  developed  fast 
prediction  algorithms  for  multiscale  graphical  models  in  spatial  and  temporal  settings.  Huang  et 
al.  (2002)  extend  the  recursive  Kalman-filter-prediction  algorithm  to  the  class  of  multiscale  tree- 
structured  models.  A  paper  giving  an  extension  of  these  spatial  tree  models  to  spatio-temporal 
tree  models  is  in  preparation.  Shen  et  al.  (2002)  use  a  multiresolution  wavelet  basis  to  detect 
changes  in  maps  from  one  time  point  to  another.  With  regard  to  statistical  theory,  Huang  and 
Cressie  (2000)  and  Cressie  et  al.  (2000)  investigate  limiting  properties  of  Markov  models  on 
directed  graphs,  which  are  called  partially  ordered  Markov  models  (POMMs).  Cressie  and  Liu 
(2001)  investigate  the  equivalence  of  certain  binary  POMMs  and  symmetric  Markov  random 
fields,  and  they  give  conditions  for  their  equivalence.  In  the  case  of  equivalence,  the  binary 
POMM  allows  fast  simulation  of  the  symmetric  Markov  random  field. 

Research  has  continued  in  exploring  the  nature  of  spatial  (and  temporal)  hierarchical 
models.  Cressie  and  Mugglin  (2000)  and  Mugglin  et  al.  (2002)  take  a  Bayesian  approach  in 
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developing  spatio-temporal  models  for  infectious-disease  count  data.  These  models  estimate 
how  an  epidemic  spreads,  spatially  and  temporally,  given  aggregated  small-area  coimt  data, 
using  the  methodology  of  Markov  chain  Monte  Carlo.  Hrafhkelsson  and  Cressie  (2003)  also 
take  a  Bayesian  approach  for  the  modeling  of  spatial  coimt  data,  with  application  to  nuclear  fall¬ 
out  data.  In  modeling  the  underlying  Gaussian  mean  process  (linked  to  the  count  data),  they 
compare  Markov  random  fields  and  geostatistical  models,  and  show  that  the  computationally 
faster  Markov  random  fields  give  very  similar  results  to  the  non-Gaussian  geostatistical  models. 
Stem  and  Cressie  (2001)  investigate  Bayesian  model  diagnostics  for  spatial  lattice  data,  an 
important  part  of  any  development  of  spatial  statistical  models.  Kaiser  et  al.  (2002)  develop 
inference  for  a  class  of  spatial  mixture  models  that  use  exponential-family  conditional 
distributions  to  model  the  underlying  latent  spatial  process. 

Impact/ Applications 

We  have  developed  two  web  sites  at  The  Ohio  State  University: 
www.stat.ohio-state.edu/  ~C2 
considers  probability  and  statistics  in  Command  and  Control. 
www.stat.ohio-state.edu/  ~sses/ research_onr.htnil 
considers  all  research  conducted  with  the  support  of  the  Office  of  Naval  Research. 
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